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Abstract 

We investigate tliree-dimensional magnetohydrodynamics turbulence in the presence of velocity 
and magnetic shear (i.e. with both a large-scale shear flow and a non-uniform magnetic field). By 
assuming a turbulence driven by an external forcing with both helical and non-helical spectra, we 
investigate the combined effect of these two shears on turbulence intensity and turbulent transport 
represented by turbulent diffusivities (turbulent viscosity, a and /? effect) in Reynolds-averaged 
equations. We show that turbulent transport (turbulent viscosity and diffusivity) is quenched by 
a strong flow shear and a strong magnetic field. For a weak flow shear, we further show that 
the magnetic shear increases the turbulence intensity while decreasing the turbulent transport. In 
the presence of a strong flow shear, the effect of the magnetic shear is found to oppose the effect 
of flow shear (which reduces turbulence due to shear stabilization) by enhancing turbulence and 
transport, thereby weakening the strong quenching by flow shear stabilization. In the case of a 
strong magnetic field (compared to flow shear), magnetic shear increases turbulence intensity and 
quenches turbulent transport. 

PACS numbers: 47.27.Jv,47.27.T-,47.65.-d 
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I. INTRODUCTION 



Most geophysical and astrophysical bodies are composed of electrically conducting fluids 
(liquid iron for the Earth and plasma for stars, interstellar medium, etc.). The evolution of 
magnetic fleld B and velocity U in these systems are often described by a simplifled model 
given by the incompressible magnetohydrodynamics (MHD) equations 

dtV + V-VV = -VP + B • VB + uAV + f , (1) 
dtB + V-VB = B ■ VU + 7]AB , (2) 
V ■ U = V • B = , (3) 

where B is the Alfven velocity (B = ISA/y/JIp where M is the magnetic fleld measured in 
Tesla; therefore B has the dimension of velocity), p is the total (hydrodynamical + magnetic) 
pressure and z/ and i] are the molecular viscosity and diffusivity, respectively. Eq. ([1]) is the 
Navier-Stokes equation including the Lorentz force which describes the effect of the magnetic 
fleld on the velocity fleld and an external forcing f , which is assumed to be at small scales. 
Eq. ([2]) describing the evolution of the magnetic fleld is called the induction equation and 
can be derived from Maxwell's equations and Ohm's law. 

In most astrophysical objects, velocity and magnetic flelds are observed to exist on a 
broad range of length and time scales. In order to characterize the evolution of flelds on 
these scales, theories, such as mean-fleld dynamo , decompose the flelds into a mean 
and fluctuating parts and parameterize the effect of the small-scale (unresolved) flelds on the 
large scale flelds in terms of transport coefficients. Specifically, expressing B = (B) + b and 
U = (U) + u, where the (•) stands for an average on the realization of the small-scale fields, 
substituting this decomposition into Eqs. ([1]13]) and averaging yield the following mean-field 
equations: 

dt{\J) + (U) ■ V(U) = -V(j9) + (B) • V(B) + i/V'(U) - V ■ 5 , (4) 
dt{B) + (U) ■ V(B) = (B) • V (U) + r]V^{B) + V xS , (5) 
V • (U) = V ■ (B) = . (6) 

The main challenge is to express the stress tensor S = (u ® u) — (b ® b) in Eq. (jlj) and the 
electromotive force £^ = (u x b) in Eq. ([5]) in terms of the large-scale variables (U) and (B). 
In the absence of non-diffusive fluxes, the stress can be expressed as Sij = —N'[j^idk{Ui) where 

2 



M^j^^i is called the turbulent viscosity tensor, which can add to the molecular viscosity in Eq. 
(jl]). On the other hand, the electromotive force is usually assumed to depend linearly on the 
mean magnetic field and only the two first term are kept (proportional to the magnetic field 
and its first derivative). Following Radler and Stepanov [3], this expansion can be written 
as: 

£i = ai,{Bj) + [/i X B]i - A,(Vx (B)), - [5x(Vx (B))], - Qjk[Vj{Bk) + Vfc(fi,)]/2 . (7) 

The first term on the right-hand side (RHS) of Eq. ([71) is the a effect which can be shown to 
generate magnetic field on large scale for a helical turbulence. It is thus a perfect candidate 
to explain magnetic fields in systems infiuenced by Coriolis force such as in the stellar 
convection zones. The second term on the RHS describes a transport of magnetic fiux by 
turbulence. The third and fourth term in Eq. ([7]) can be described by introducing an 
anisotropic turbulent diffusivity. The last term proportional to n does not allow a simple 
interpretation. The presence of the additional terms besides the a and /3 effect is possible 
only for anisotropic or/and inhomogeneous turbulence. 

There has been accumulating evidence that a strong shear reduces turbulent transport via 
shear stabilization by flow shear |4|]. This is basically because shear advects turbulent eddies 
differentially, elongating and distorting their shapes, thereby rapidly generating small scales 
which are ultimately disrupted by molecular dissipation on small scales (see Fig. [I]). That 
is, flow shear facilitates the cascade of various quantities such as energy to small scales (i.e. 
direct cascade) in the system, enhancing their dissipation rate. As a result, turbulence level 
as well as turbulent transport of these quantities can be signiflcantly reduced compared to 
the case without shear. Another important consequence of shearing is to induce anisotropic 
transport and turbulent level since flow shear directly influences the component parallel 
to itself (i.e. the x component in Fig. 1) via elongation while only indirectly the other 
two components (i.e. the y and z components in Fig. 1) through enhanced dissipation 
and incompressibility. Indeed, the flow shear has been shown to signiflcantly reduce the 
turbulence intensity and the turbulent transport of angular momentum, particle mixing and 
magnetic diffusion both in hydrodynamics |5|, l6| and magnetohydrodynamics 7]]. In 81], by 
assuming a uniform large-scale magnetic fleld, we showed that the a effect is quenched by 
flow shear and magnetic fleld. 

In this paper, we flrst examine the effect of non-uniform magnetic fleld on the linear 
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FIG. 1: Sketch of the effect of shear on a turbulent eddy. 



stability of MHD fluids. We tlien proceed to investigate tlie effect of a non-uniform magnetic 
field on MHD turbulence. A consistent derivation of turbulent magnetic diffusion [jS effect) 
of non-uniform field in 3D MHD is crucial since there has been controversy over the f3 
suppress,™, by strong .nagnetic field fi. In companson. in 2D MHH the effect was shown 
to be severely quenched by magnetic field and also by shear flow [10|. Furthermore, in the 
presence of shear flow, another interesting effect of magnetic shear in non-uniform magnetic 
field is that the latter can interfere with flow shear and thus weaken the quenching of 
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3. 



We investigate 



turbulence by flow shear, as shown in 3D reduced MHD turbulence 
if this result holds in 3D MHD by quantifying the effect of magnetic shear on turbulence 
intensity and turbulent transport. The reminder of the paper is organized as follows. In 
Section [ni we present the main governing equations and the quasi-linear approximation that 
is used to solve for the turbulent fields. We then present the calculations of the transport 
coefficients (turbulent viscosity and a effect) in the limit where the uniform component of 
the large-scale magnetic field is negligible in Sec. IIIII and then in the other extreme limit 
where it is very strong (compared to shear flow) in Sec. IIVI Conclusions are provided in 
Sec. lYl 



4 



II. MODEL 



A. Magnetohydrodynamical equations 

We use the similarity of the Navier-Stokes ([T]) and induction equation ([2]) and work in 
terms of Elsasser variables 13j instead of the (fluctuating) velocity and magnetic fields: 
= U + B and = U — B. Assuming a unit magnetic Prandtl number {u = i]), the 
MHD equations ([T]l3]) can be rewritten as: 

9t*+ + ^ ■ V*+ = -VP + z/A*+ + f , (8) 
dt"^- + ■ V*" = -VP + uA<b + f , (9) 
V ■ = V ■ = . (10) 

To study the effect of shear flows and magnetic fields on small-scale turbulence, we prescribe 
a large scale flow of the form (U) = —xAsy and a non-uniform large-scale magnetic field 
(B) = {Bq — Bix) By. The chosen configuration (with parallel velocity and magnetic field) 
ensures that there is no direct influence of (U) on (B). To solve the equations for the Elsasser 
variables, -0+ = - = - (U + B) and ^ = <Sf - = - (U - B), we 

use the quasi-linear approximation assuming that the interaction between fluctuating fields 
is negligible compared to the interaction between large and small-scale fields. The equations 
for the fluctuating fields can then be written as: 

+ VV^ + V"- V(*+) = - Vp + z/AV'+ + / , (11) 

dti/^- + {<if+) -Vjp- + ip+ -Vi^-) = -Vp + uAil^- + f , (12) 

V • = V ■ = , (13) 

where p is the fluctuation in the pressure. To solve these equations, we use vanishing initial 
conditions: i/^^ = i/'^ = at t = to- An equilibrium is then reached at long times when the 
power injected by the forcing balances the dissipation. To solve Eqs. flllH13l) . we introduce 
a time-dependent Fourier transform: 

y(x,t) = j ci3kg4*^-W-+'=-f+'=^#(k,t) , (14) 

where we choose k^it) = kx(to) + Aky(t — to) to account non-perturbatively for the effect of 
the non-uniform components of (U) and (B). Fourier-transformation of Eqs. ()1 111131) leads 



to the following equations: 



dti^f T kyBidkJt -{A± Bi) i,f5i2 = ±iB^kyi,f - ihp - uk^f + h , (15) 

+ kyi^^ + = . (16) 



Note that the time-dependence of the wave number cancels exactly the advection by the 
mean velocity shear (see Kim for details) and that the second term on the left hand side 
of Eq. (ITSl) is obtained by Fourier transform of the advection by the mean magnetic shear 
in the following way: 

rr{B^xdyi}f) = j d^xe-'\-''^^'^''~"'«y^^^'\Bixdy^t (17) 



= -kyBidkM ■ 

Changing the time variable from t to r = k.j.{t)/ky = kx{to)/ky + A{t — to), the first two 
terms on the left-hand side of Eq. (IT^ can be grouped together as dt = Adr and kydk^ = dr- 
Assuming the forcing to be incompressible, Eqs. ( |T5lfT6l) can be rewritten: 



;i T n)dr^f - (1 ± n)^j5i2 = 



(18) 



g2 _|_ ^2 



;l-7^)^+ + (l + 7^)^^ 



1 ~„ 







(19) 



Here, IZ = Bi/ Ais the ratio of the magnetic shear to the velocity shear; 7 = B^ky/A is the 
ratio of the Alfven frequency to the fiow shear; b = k^/ky and (7^ = 1 + 6^; ^ = vk"^/ A. Note 
that Eqs. (11811191) are invariant under the following transformation: 7 —7, TZ ^ —TZ and 
^ i/j^. Consequently, can be obtained from ij)^ (and vice versus) by changing the 
sign of 7 and TZ. Using the following variables: 



the x-component of Eq. ( |T8l) can be rewritten: 

^70 



and iIj^ (r) 



(20) 



1 — it Q'^ + T'^ 



(21) 



The coupled equations for 0+ and in Eq. (12T|) can be combined to a closed equation for 



dl 



+ 



+ 



T 1 - 7^2 

il [ 1 



r 



f22) 



(^2 + ^2) V(i -7^)r l + 7^y l-7^ 



^2 27^r 
— + - 



7^ 



+ 
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l-7^2 



^7 + ^{g^ + 

1 + 7^ 



Two initial conditions are needed in order to solve Eq. fl22l) . At r = tq = kx{tQ)/ky, we 
assume that initially there is no velocity and magnetic perturbations (0^(t = tq) = (/)^(r = 



^0 



0). The use of < 



^0 



in Eq. (EI]) gives us [9,0+] (r = tq) = /.(r = ro)((7' + r2)M 



as the second initial condition. As Eq. f l22i) cannot be solved in general case, we will 
consider the two cases of weak and strong magnetic field, given by 7 = B^ky/A ^ 1 and 
7 3> 1 respectively. Note that the weak magnetic field limit does not restrict the magnitude 
of the magnetic shear. The only constraint on the magnetic shear comes from our assumption 
that the system is stable (i.e. \TZ\ < 1 as shown in the next subsection). 



B. linear Stability 



We start by studying the stability of the large scale fields (U) and (B) by considering 
the behavior of the perturbations of Eq. (1221) in the long time limit. As Eq. (I22p is a second 
order differential equation in r, the homogeneous equation has two independent solutions 
01 (r) and 02 (t). Using WKB theory (see Appendix [Bl for details), we can show that these 
two functions have the following asymptotic behavior in the large r limit: 



exp 



t>2[T) ~ rexp 



1 + 7^ 



Qir) 



Z7 



Qir) + 



l + 7^ 

27 



-T 



(23) 
(24) 



where Q{x) = g'^x + x^/3 and the ~ symbol stands for asymptotic behavior in the large r 
limit. Eqs. (12311241) show that if \TZ\ > 1 one of the solutions of the homogeneous problem 
is exponentially divergent for large r. This is a similar result as that found by 1^ who 
studied the resistive tearing instability in the context of fusion plasmas. Interestingly, the 
tearing instability is stabilized by flow shear if \TZ\ < 1. Note that the limit 7^ is 
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singular as we do not recover the result of [5|] with no magnetic shear. This is due to an 
additional dissipative layer imposed by the magnetic field. Consequently, the effect of the 
magnetic shear, is two-fold. On the one hand, it renormalizes the diffusion rate ^ in the 
exponential factor to .^/(l — 7?.) and ^/{l + TZ) (with a continuous limit to the case TZ = 0). 
On the other hand, magnetic shear changes the scaling with r in front of the exponential 
in a non-continuous fashion (note that the limit ^ = is also singular as this asymptotic 
behavior does not recover the result for the ideal system that can be computed exactly). In 
the following, we restrict our study to the case where the system is stable {\1Z\ < 1). 



C. Transport coefficients 

Our main interest is in the total stress and the electromotive force, which determine the 
growth/decay of the large-scale velocity field and the large-scale magnetic field, respectively. 
The assumption of a large scale fiow is in the y direction and depending only on the x 
coordinate has two implications for the stress. First, only the components Mfj^y do not 
vanish. Second, only the y-component of Eq. (jlj) is of interest and the divergence of the 
Reynolds stress reduces to dxSxy (the two other terms involves derivative with respect to y 
and z). Consequently, in the following, we are interested in only one component of the total 
stress S defined as: 

S = S^y = {U^Uy) - {hj)y) = ^ ( ^ ^ + ^ ) " (25) 

Note that this total stress consists of the difference between the Reynolds stress (u^Uy) 
and Maxwell stress {hj)y). In the following, we refer to turbulent viscosity as the only 
component of interest v'^ = Af^yxy For the assumed shear flow (U) = —Axey, the turbulent 
diffusivity can be computed as S = v'^ A. 

Similarly, as we chose the large-scale magnetic fleld to depend only on x, the only com- 
ponents of V X £^ are the one in the y and z direction given by: 

Sy = {u.K - u^K) = - CV'+)/2 , (26) 

= {uxby - UyK) = {^P^^p- - ry^t)n ■ (27) 

For our chosen configuration of the magnetic field which depends only on x, the electromotive 
force has only the following terms (see jisl for a general expression of the electromotive force 
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for arbitrary shear flows): 



S-y = OiyyBo , (28) 

= {a,y + /i,)5o + , (29) 

where ayy and a^y are the components of anisotropic (due to shear flow) a effect, /i^ is the 
turbulent transport of magnetic flux and (3 characterizes the turbulent (3 effect. Note here 
that only these three coefficients are present in our configuration. In particular, phenomena 
such as the $7 x J 3| and shear current effects [l^, which have been advocated for turbulence 
affected by rotation and shear, are absent here. For instance the shear current effect vanishes 
because the large-scale vorticity and the curl of the magnetic field are parallel to each other 



D. Forcing 

To calculate the correlation functions involved in the transport coefficients [see Eq. ( 1251 - 
[271) ]. we consider an incompressible forcing which is spatially homogeneous and temporally 
stationary with a short correlation time r/. Specifically, in Fourier space, the correlation 
function of the forcing is taken as: 

{/^(ki,ti)/,(k2,t2)) = Tf (27r)35(ki + ka) 6{h - t2)f^^,{k^) . (30) 

As noted previously, the a effect can be linked to the helicity of the turbulent flow. Conse- 
quently, we consider a forcing with both a symmetric part (with energy spectrum E) and a 
helical part (with helicity spectrum H) given by: 



.(k) = Eik) - ^) + iei^,k,^ . (31) 



In the following, the turbulence intensity, turbulent viscosity and a effect are expressed in 
terms of the kinetic energy cq = (u^) and helicity = (u(-V x u) of the flow created by 
the forcing / in absence of shear and magnetic fleld: 

-i^r*^' (3.) 
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The presence of the viscosity v in the denominator of Eqs. fl32H33p is due to the fact that 
without dissipation (z/ = 0), there is nothing to dissipate the energy injected by the (small- 
scale) forcing causing the accumulation of small-scale fields. Therefore, the growth of kinetic 
energy or transport is unbounded. In appendix [XJ we show the derivation of Eq. fl5^ and 
Eq. dMD. 



III. WEAK MAGNETIC FIELD (7 < 1) 

In this section, we investigate the influence of magnetic and flow shear on turbulence 
properties in the limit of a very weak magnetic field (i.e. 7 <^ 1). In the following, we 
thus keep the magnetic field only to the lowest order to compute the a effect. As it seems 
impossible to exactly solve Eq. fl2^ even for 7 = 0, we obtain our results in two different 
limits of weak = uky/A ^ 1) and strong <^ 1) shear in comparison with the diffusion 
rate (i^ky). 



A. Weak Flow Shear: ^ > 1 



In the case of a shearing rate much weaker than diffusion rate = uky/A ^ 1), we can 
perform a WKB analysis (very similar to that of appendix [B| of Eq. (l22l) by using ,^ ^ 1 as 
a large parameter, obtaining the following two solutions of the homogeneous system ( l22l) : 



(^2 + r2)3/2 



exp 



'^2[T) = A/fi'^ + t"^ exp 
to leading order in Here: 



1 + 7^ 



Qir) 



«7 



-T 



1 + n 



n 



- (g2 [E- 

\J g"^ + t"^ exp \E~^{t 



E~{x) = -e/(l + 7^)g(x) -27/(1 + 7^)a; , 

E'-ix) = -e/(l-7^)g(x) + 27/(1 -7^)x. 



(34) 
(35) 

(36) 
(37) 
(38) 



Using the method of variations of parameters, the full solution of Eq. (1221) with the initial 
conditions 0"'"(r = tq) = and [dr(t>'^]{T = tq) = /^.(r = ro){g'^ + Tq)/A is found as: 

{9' + n')Uu)^/^^T^ 



du- 



A 
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exp [£' 



E+{u 



(39) 



where is defined in Eq. ( 138|) . The z component of the field can be obtained by integrating 
the z component of Eq. f|T8l) with the following result: 

-6£rf^ V^^|[^| Mexp - /+(^,r) . (40) 

Here we defined the integral /"*" as follows: 

(41) 



r dx 



:{Q(a:) - Q(^)} - T^^i^ - u) 



l-7^2'^' ' l-7^2 

The y-component can be obtained using the incompressibility condition. 

Using Eqs. ( l39M0l) . we can compute the turbulent intensity and turbulent transport 
coefficients in a similar way as shown in the beginning of Appendix [El The results are as 
follows: 

^^'^ = (2^ / d'k'^E{k)U^) , (42) 

(bl) = / d'k'^Eik)lUk) , (43) 

(nl) = I d'kE{k)lUk) , (44) 

(bl) = 777^ I d'kEik)h,{k) , (45) 



(27r)3^ 

where kjj = ky + k1 and the various integrals J's are provided in Eq. (ICltlCTp of Appendix 
[Cl In the limit ,^ ^ 1, the approximate value of these integrals in Eqs. fl42ti48p can be 
evaluated and then used for the computation of the turbulent intensity and transport. The 
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results are: 



^ ^' (27r)3 y 2uk^ 1 - 7^2 ' ^ 
^^^^ ~ (2vr)3 J ^ " 2vk^ 1 - 7^2 ' ^^"^ 
^""^^ ~ (27r)3 j ^ 2uk^ 1 - 7^2 ' ^^^^ 



(27r)3 y 2z/A;4 1 - 7^2 ' 



Using the fact that the forcing is isotropic, Eqs. (1491155^ can be simphfied by integration 
over the angular variable 6 and (after expressing the wave vector in spherical coordinates: 
kx = kcos6, ky = /c sin 6' cos and kz = fcsin^sin0). Eqs. (14911551) can then be recast as: 



(bl) = 7?^, rdk^-^^^^ ~ eoW,(7^) , (59) 



(27r)2 7o 3z/ 1 - 7^2 

r/ ,,H{k) ho 

^^(£,r-^^^^^-^^^ 

In Eqs. (I56ti62p . the functions W„, which are plotted on Fig. [2], characterize the dependence 
of turbulence intensity and transport on the magnetic shear. Figure [2] shows that the 
magnetic shear tends to increase the turbulence intensity whereas it decreases the turbulent 
dissipation of both momentum and magnetic field {vt and (3). Furthermore, the a effect is 
not affected by the magnetic shear and is the same as in the kinematic regime. The increase 
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FIG. 2: Dependence of yVi{R), W2(i?) and W^iR) on the magnetic shear. 

in the turbulence intensity can be understood in terms of instability of the homogeneous 
part of Eq. (!22l) . As can be seen in Eqs. fl23H24l) . the system tends to becomes unstable 
(the exponential decay being slower) when the magnetic shear is increased. These results 
are summarized in Table [B 




B. Strong Flow Shear: ^ < 1 



The homogeneous part of Eq. fl22|) has two independent solutions 0i(t) and 02 (t). 
Guided by the result of section IllBt we rewrite these two solutions as: 



01 (r) = Ci(r)exp 

02 (r) = C2(r)exp 



1 + 7^ 



Qir) 
Qir) 



«7 



-r 



1 + 7^ 

i'-f 



(63) 
(64) 



i-n^' ' l-7^ 

We can determine the asymptotic behavior of the unknown functions: Ci(r) ~ and 
C^ij^ ~ T and for r ^ 1 (see appendix [B] for details). Using the method of variation of 
parameters, the general solution of Eq. ( !22l) is obtained as: 



du 



TO 



A 



u 



u "" \ W{u) 



C2(r)Ci(M)e^^(^)+^"(") - C2(ix)Ci(r)e^+(")+^'(^) j 



il + i{9^ + m') C2(T)Ci(u)e^+W+^"(") - C2(M)Ci(r)e^+(")+^"W 



1 + 7^ 



W{u) 
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Here W{t) is the Wronskian of the two solutions of the homogeneous problem; £"+ and E~ 
are defined in Eqs. (13711381) . After a long algebra, we can compute all the components of the 
fields to leading order in ^ with the following results: 

0+(r) ~ r du ^^'^'^y''^''^ b3NC2(r)e^^(-)-^^(") + C4(^)Ci(r)e^-(-)-^-W] (66) 

^+(r) ~ r du-J^^e^"-^-^-^^^-^ 



r du ^^ +^ U+M-i^+W + ei^-M-i^-Hl c,{u, r) . (67) 



i-n 

Here, note that the analytical form of the functions C^{u)^ C^{u) and C^{u, r) are unknown. 
However, their asymptotic behavior for large arguments (r oo) can be computed from 
those of Ci and C2. In particular, in the large r limit, the functions C5 is independent of 
T. As shown in appendix [El the asymptotic dependence on r is sufficient for computing 
the scaling of the turbulence intensity with the shear. By using similar techniques as in 
appendix [El the magnitude of the turbulent velocity and magnetic field can be found to 
leading order in ^ as: 



i<) = .0 w/' / rf^fc^-^f^J.(k) ~ eo^V.m , (68) 



(27r)3(l - 7^2)M 7 k 



{<) = ^^^^ I d^kEm-'" [^.l(k)V2(7^) + J.2(k)V3(7^) + J.3(k)V4(7^)] (69) 

~ eoe'/'[V2(7^) + V3(7^) + V4(7^)] , 
{^D = .o„^.v7^L.^. . / d'k'^^^JM eoiUm , (70) 



(27r)3(l - 7^2)M 7 

{hi) = ^^^^Jd'kE{k)r'^'[J,,{k)M2{n) + J.^^^^ (71) 
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Here, the J's are convergent integrals which are independent of the velocity; V's and A4's 
are defined as: 



Vl(7^) 

V2(7^) 
V3(7^) 
V4(7^) 

Miin) 
M2{n) 



1 



A [(1 _ 7^)-5/3 + (1 + 7^)-5/3 ^ 2(1 - 7^2)-2/3] , 

i(V2(7^) + (l-7^2)-^/^) , 



[(1 - 7^)l/3 + (1 + 7^) V3 + 2(1 - 7^2)l/3] 



4(1 -7^2)2 

7^2 
(1 -7^2)2 ' 

i[(l-7^r^/3 



(l+7^)-5/3_2(l_7^2)-2/3] ^ 



-(7W2(7^) + 7^2(l-7^2)-'/') 



[(1 - 7^)l/3 + (1 + 7^)l/3 - 2(1 - 7^')^/^] 



(72) 

(73) 
(74) 
(75) 

(76) 

(77) 
(78) 
(79) 



4(1 -7^2)2 

Eqs. ( I681I7T]) show that both the turbulent intensity and the magnetic field are reduced by 
strong fiow shear A. Furthermore, the quenching is anisotropic as the components in the 
direction of the shear ((f^) ~ (&^) ~ A~^) are much more reduced than the components 
in the perpendicular direction {{vD ~ {hi) ~ A""^^^) as A increases. Figure [3] shows the 
dependence on the magnetic shear TZ of the velocity and magnetic field amplitude: they 
are increasing functions of the magnetic shear. This shows that the effect of the shear is 
to increase the turbulent intensity. In other words, the magnetic shear acts in the opposite 
way to the velocity shear, interfering with fiow shear to weaken the quenching of turbulence 
by fiow shear. 

In the limit ^ <^ 1, the Reynolds stress and the electromotive force can similarly be 
computed with the following turbulent transport coefficients: 



T 



/5 



(27r)3^2 
(27r)3^2 



d^kE{k) 



l-7^2 



d^kE{k)I,3(k)- 



7^2 



eeo 



;i -7^2)z/A;2 



In 



7^2 7^ 



In 



1 + 7^ 



eoe 



l-7^2 
1 



;i -7^2)^/p7^ 



In 



7^ 



7^ 



(80) 



(81) 



Eq. (]80l) shows that the turbulent viscosity and diffusivity are reduced by a strong shear. 
Moreover, the reduction in the turbulent viscosity is weaker by a logarithmic factor than 
that in the turbulent diffusivity. Eq. fl63l) also shows that the effect of the magnetic shear 
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FIG. 3: Functions characterizing the dependence on the magnetic shear of the velocity (left) and 
magnetic field (right). 

is to increase the turbulent viscosity and diffusivity. It thus weakens the quenching effect of 
the shear flow. These results are again summarized in Table ([T]). 



IV. STRONG MAGNETIC FIELD 



[ ^/^^T^^{u)exp[E+{r)-E+{u)]du, (82) 

•J TO 



For strong magnetic field (7 = B^ky/ A ^ 1), a WKB analysis of Eq. (I22l) gives: 
^ 

^(l-7^)VF+^ 

I exp[E+(r)-E+(n)]/(n,r)V7T^/,(n)dn, (83) 

J TO 



h 



to leading order in 7"^. Here: 

du 1 I T u 



f 

J u 



(84) 



and and E~ are given in Eqs. (13711381) . 

The solution for the conjugate Elsasser variables, and ^ can be obtained by changing 
the sign of 7 and TZ in Eqs. ( I82ll83|l . Using Eqs. ( I82ll83l) . we can compute the intensity of 
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turbulence, with the following result to leading order in 7 ^: 

^"^^ = ^^^^ = W^lWi^ J ^T^ ^ ^ (T^W ° ' ^ ^ 

1 /3 

{<) = m = I d'kE{k)S'{Rni/3) (I) ~ e^'S\R)eo . (86) 

Here, kjj = ky + k^ and T(k) = Ti\ky\/2kH — arctan(A;^//i;//). The fact that the velocity 
and magnetic fields are in equipartition is due to a strong magnetic field which drives an 
Alfvenic turbulence. Note that the effect of the magnetic shear is included in the term: 

S%R) = ^ [(1 - 7^)-^/=^ + (1 + 7^)-^/=^] . (87) 

Eqs. ( I85ll86|) show that the turbulence is unaffected by Bq while severely quenched by flow 
shear. Furthermore, turbulence is less severely reduced in the direction perpendicular to the 
shear than in the direction of the shear. In particular, we see that, even in the limit of strong 
magnetic field, the turbulence intensity is reduced solely by the flow shear. This is because 
the magnetic field forces the turbulence to be more wave-like (due to Alfven waves), thus 
increasing the memory time, without necessarily reducing its amplitude. Note that both 
(m^) and {uD increase with the magnetic shear. That is, the effect of the magnetic shear is 
again to weaken the quenching of the flow shear. 

Similarly, the turbulent viscosity and the a effect are found as: 



" 2(27r)3 J """^ k^Blkl + z/2A;2) ' 
to leading order in 7"^. The a-effect is the same as in the case without magnetic shear 
showing that the magnetic shear has no effect on the a-effect. In comparison, the turbulent 
viscosity z/^ is reduced as the magnetic shear increases. Furthermore, we find that the 
turbulent diffusivity (3 vanishes to leading order. Non-trivial diffusivity will be found only 



at the higher order in 7 ^. This means that magnetic fields 
too strong, which is in agreement with numerical simulations 



lardly diffuse when they are 
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V. CONCLUSION 



To understand the properties of astrophysical and geophysical magnetic fields, we inves- 
tigated the turbulent transport of non-uniform magnetic field in the presence of flow shear. 
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Weak flow shear ^ = uky/A » 1 


Strong flow shear ^ 1 
































(3 







TABLE I: Summary of our results obtained in the weak magnetic field limit. f+{R) and f~{R) 
represent a function which increases or decreases with magnetic shear TZ = Bi/A, respectively. 



The linear analysis reveals that magnetic shear stronger than flow shear leads to instability. 
We considered the two limits of strong and weak magnetic field. In the weak magnetic field 
limit, the magnitude of the magnetic shear is limited only by the stability condition that 
the magnetic shear is to be weaker than the flow shear. 

In the case of a weak magnetic field, the flow shear is shown to reduce both the turbulence 
intensity and turbulent transport due to shear stabilization in agreement with our previ- 
ous results [5], ^. In particular, turbulent viscosity and magnetic diffusivity {(3 effect) are 
strongly suppressed as A'"^ for strong flow shear A. When magnetic shear is incorporated 
the turbulence intensity is increased. This is due to the fact that magnetic shear tends to 
make the system unstable as shown by our stability analysis. In contrast, the turbulent 
transport (both turbulent viscosity and diffusivity) is reduced by magnetic shear if the flow 
shear is weaker than diffusion rate while it is increased if the shear is stronger. In all the 
cases, we found that, for strong shear, the magnetic shear opposes the effect of flow shear, 
thereby compensating the quenching of turbulence by shear stabilization. These results are 
summarized in Table [B 
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For strong magnetic field, we found tliat tlie turbulent intensity is not quenched by strong 
magnetic field. In comparison, the turbulent viscosity and a effect are reduced by strong 
magnetic field with scalings B^^ and -B^^, respectively. Furthermore, we found that the 
turbulent diffusivity of magnetic field (the (5 effect) vanishes for strong magnetic field to 
leading order, indicating that strong magnetic fields hardly diffuse. To recapitulate, the /3 
effect is quenched by magnetic field for a large constant magnetic field whereas in the weak 
magnetic field limit, it can be reduced by strong flow shear. 

It will be interesting to extend our theory to incorporate the effects of rotation which 
will consistently give rise to a effect and non-diffusive momentum transport (A effect) due 
to shear-induced anisotropy How the A effect, a effect, turbulent viscosity and particle 
transport are affected by rotation, magnetic field and shear would be of great interest with 
important practical implications. These issues will be addressed in future publications. 
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APPENDIX A: KINETIC ENERGY AND a EFFECT IN THE ABSENCE OF 
SHEAR 

In the case without large-scale shear flow, the linearized equation for the fluctuating 
velocity can be written as: 



In the case where the forcing is incompressible, the pressure vanishes (p = 0) and the solution 
of this equation can easily be obtained in Fourier space as: 




(Al) 




(A2) 



19 



Using Eq. (lA2p . the turbulent intensity can be computed as follows: 



(27r)^ J Jo Jo 



(A3) 



(27r)3 j ' ukl 



{2n) 



3 



The second line in Eq. ( \A3\} is obtained by using the correlation function of the forcing 
given by Eq. (!30ll . The fourth line is obtained by taking the long-time limit {t oo) while 
the last line is obtained by integrating over the angular variables. 
Similarly the kinetic helicity can be obtained as: 

(u-(Vxu)) = — ^ / d^kid^k2e'^^'+^^>'' [ dui [ rfM2 e-'^^'=?(*-"^)+'=2(t-2)} X 
(271") J Jo Jo 

{eimp {ik2m)fl{^l,Ul)fp(i^2,U2)) (A4) 

= j d^h ^ dm e-2-^''(*-«^)Q™p(-zfci™)Kzp(-ki) 

= f d^ki f dui e-^-'^?^*-"^) 2H{ki) 

(27r) J Jo 



(27r)3 J ' ukl 

2-/ r°^.. H{k 



(2vr) 



Jo V 



APPENDIX B: WKB ANALYSIS IN THE LONG TIME LIMIT 

To study the behavior of Eq. (1221) for large time, we introduce a small parameter e and 
write T = x/t. In terms of the new variable x, the homogeneous part of Eq. fl22|) can be 
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rewritten as: 



e 2^ 



X 



9 + ^ 



+ 



«7e 



76 



(^2e2 + 3,2) _ 7^)2; e(l + 7^)y l-7^V ea; e(l + 7^) 



27^x 



+ 



l-7^2 



Using the following WKB ansatz: 

(f)~^{x) = exp 



{So{x) + eS^ix) + 



(B2) 



and solving order by order (in e), we can obtain the two solutions of the homogeneous 
equation (IBip with the following values of S'o, Si, S2 and S3: 





Solution 1 


Solution 2 


So{x) 






3(l+7^) 


3{l-1l) 


Si{x) 








S2{X) 






S^ix) 


—2 Inx 


Inx 



Plugging this result into Eq. (\B2\i and changing back to the original variable r = x/e, we 
obtain the two following approximate solutions to the homogeneous equation: 

. , . 1 



exp 



62 (r) ~ rexp 



1 +7^ 



1 + 7^ r 



-Q(r) + ^^r+ -/2(r) 



(B3) 
(B4) 



Here, Q{x) = g^x + h and Z2 are two functions which converge in the large r limit. 
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APPENDIX C: INTEGRALS IN THE WEAK SHEAR LIMIT 



In the computation of the turbulent intensity and transport coefficients, we obtained the 
following integrals in Eq. 



+00 



dr 



i-bx 



A- -I- pos 

, 4((72 + r2) I (1 -7^)2 ^ (l + 7^)2 ^ l-7^2 

+00 j g2_B+ (a,T) ^E^ia^r) 2e'^^0^'^''^^ 

+ : ::z7r- COS 



27 



„ 4(^2 + 7-2) (1 _ 7^)2 (l+7^)2 1_7^2 



1+^2 ( ^Eoia,T) (,2E+(a,T) 2e2'^o("'^) 



+ 



+ 



4(^2 + a2) (l-7^)2 ' (l+7^)2 ' l-7^2 



■ COS 



l-7^2 

27 

l-7^2 

27 



(r-a) 



r - a 



l-7^2 



(r-a) 



ai,2 [ gi?o+(a,T) g2i<;o-(a,r) Q2Eo(a,T) 



y/^T^ [(l-7^)2 ^ (l + 7^)2 ^ ' l-7^2 

52 „2 r ^E+(a,T) ^ p2£;o-(a,T) 



j^l fe^^^(-~-)/(^-'^')[I+(a,r) + r{a,r)] 



o2Eo{a,T) 



^2 + a2 (1 - TZy 



(l + TZ) 



+~ l + a^ r eEo-(a,T) (,2E+(a,T) 2e'^^°^'''^^ 



dr 



+ 



4(^2 + ^2) (l-7^)2 (l + 7^)2 l-7^2 



COS 



27 



1 -7^2 



(r-a) 



52 [ gSg+Ca.r) ^2Eg (a,T) g2Eo(a,T) 



/7 - 



y^2^ ^^(l_7^)2 ^ '(l + 7^)2^ l-7^2 

52^2 r gi?o+(a,r) 



Rl fe2^^("-")/(i^^')[/+(a,r) + /-(a,r)] 



e2So~(«.^) 



2Eo(a,T) 



g^ + a^ I (1 -7^)2 
dr 



+ 7r^l^'(«'^)l' - T-R^^^ (e^^-(— )/(^-^^)[/+(a,r)/-( 



(1 + ^) 



5(2 + r2 



2Eo{a,T) 



COS 



27 



l-7^2 



(r-a) 



r(5(2 + a^) _^ ab"^ 



(72 + r2 ^/ g"^ + cl"^ 



b\g' + a^] 



4 



'^dr-^ 



,2Eo{a,T)f^ 



sm 



r^ 



27 



1 -7^2 



.r - a 



r> + 00 



la + 
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Here, Rl stands for the real part; /"'" is defined in Eq. (HTi) and / is obtained from J"*" by 
changing the signs of IZ and 7; we defined a = kx/ky and the following functions: 



a, r 



Eq (a, T 



Eo{a,T 



lQ{a,T 



i-n 
1 + 7^ 
1 -7^2 



- g(a)] , 

[Q{r) - Q{a)] , 
[Q{r)-Q{a)] , 
1 + exp 



dx 



(^2 + 3.2)3/2 

(^2 + a;2)3/2 

(5f2 + x2)3/2 



1 + exp 
1 + exp 



1 -7^2 

1 -7^2 
1 -7^2 



{Q(x)-g(a)} 
{Q{x)-Qia)} 
{Q(x)-Qia)} 



cos 



sm 



27 



l-7^2 

27 
l-7^2 



Uc — a 



X — a) 



(C8) 
(C9) 
(CIO) 
(Cll) 

1?) 

r(Si3) 



The integrals Jq (a, r), (a, r) and 1^ (a, r) are obtained by changing the signs of 71 and 7 
in Eqs. flC8tlC13D . 



APPENDIX D: ASYMPTOTIC EXPANSION 

Expanding all the integrals in Eq. (ICip in powers of we obtain the following two 
leading orders (by keeping only the terms which are even in all wave-numbers as the terms 
with an odd number of wave-numbers would vanish upon angular integration): 

1 / I 1 \ _ 1 1 -7^2/2 

~ 8e((?2 + a2)2 [jr^ ^ (TT^ ^ ) ~ 2e(^2:f^T^^ ' ^ ' 

8e(^?2 + a2)2V(l-7^) (l + 7^) ; 2^((72 + a2)2 1 - 7^2 ' ^ ' 

8e((72 + a2)2V(l-7^) (l + 7^) ; 2e((72 + a2)2 1 - 7^2 ' ^ ' 
^ l + g^ / 1 1 _ \ ^ 1 + a2 7^2/2 

8^(^2 + ^2)2 V(l-7^) (l + 7^) j 2e((72 + a2)2 1 - 7^2 ' ^ ^ 

" 2^(^2 + ^2)2' 
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APPENDIX E: INTEGRALS IN THE STRONG FLOW SHEAR LIMIT 



As an example of how to compute the correlation functions, we show the main steps to 
obtain (X^). The Fourier transform X of X is given by: 



X = 

where Q is defined in Eq. 
[see y, for details]: 



TO 



du G{u)F{r) exp K {Qir) - Q{u)}] fi{u) , (El) 
The correlation of the variable X can then be computed as 



d^k Kii(k)G{a) 



F\T)exp[-aQ{r)-Q{a)}] dr 



(E2) 



{2tt)^A 
where a = k^/ky. 

For the computation of the correlation in the strong shear limit ^ 1), we have to 
compute r-integrals of the form: 



r+oo 

K = / F\t) exp [-e {Q{r) - Q{a)}] dr . 

J a 



(E3) 



Here F"^ has the scaling F'^{t) ~ in the large r limit. When 6 < —1, K exists for = 
thus, the ^ ^ 1 limit can be easily obtained by putting ^ = in Eq. (]E3p . This is however 
not the case when ^ > — 1 as the integral diverges as ^ ^ 0. In that case, by making the 
substitution y = 2^r'^/3, the integral ( 1E3I) can be computed in the ^ <^ 1 limit as: 



K 



+ 00 
+ 00 



2^ 



1/3 



-y 



2 \ 1/3' 



2^ 



2^ 



2^ 



(0+l)/3 



(e+i)/3 



exp 

dy 

oo 

eM-y]y^'-'^''dy 
e + i 



2^ 



2e 



(E4) 



where F is the Gamma function. In summary, if the integrand has a scaling (excluding 
the exponential factor) in the large r limit, the correlation function scales as in 
the strong shear limit (^ <^ 1). For instance, Eq. ( 1671) shows that in the large r limit, the 
integrand does not depend on r (0 = 0). Therefore, for the computation of (m^), the r 
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integral has the scahng ^ with the shear. 
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